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ABSTRACT 


Examination  of  the  two  electron  density  matrix  in  the 
alternant  molecular  orbital  approximation  for  a  large  system 
in  a  singlet  state,  reveals  that  there  are  both  short  and  long 
range  correlations  between  the  electrons  of  opposite  spin.  There 
is  also  a  long  range  correlation  between  electrons  of  the  same 
spin.  These  correlations  are  not  found  when  the  wave  func¬ 
tion  for  the  system  is  approximated  by  a  single  determinant 
of  alternant  molecular  orbitals. 
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1 .  INTRODUCTION 

One  of  the  more  practical  methods  that  has  been  suggested  for 
including  some  correlation  between  electrons  in  a  system  of  many 
electrons,  is  Lowdin's  alternant  molecular  orbital  (AMO)  method.  *  In 
this  paper,  we  show  that  using  the  AMO  method  for  a  large  system  in 
a  singlet  state,  it  is  possible  to  give  general  and  easily  interpreted 
formulas  for  the  spin  components  of  the  one  and  two  electron  density 
matrices.  Examination  of  the  spin  components  of  the  two  electron 
density  matrix,  reveals  both  short  and  long  range  correlations  between 
electrons  of  opposite  spin,  and  a  long  range  correlation  between  elec¬ 
trons  of  like  spin  in  addition  to  the  short  range  statistical  correlation, 
or  Fermi  hole.  Our  discussion  serves  to  emphasize  the  differences 

between  the  spin  properties  of  approximate  wave  functions  which  are 
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eigenfunctions  of  the  total  spin  angular  momentum  operator  S  and 
those  which  are  not. 

In  the  next  section  we  briefly  describe  the  properties  of  AMO’s 
and  then  write  out  formulas  for  the  spinless  one  and  two  electron 
density  matrices  of  a  system  of  many  electrons  in  a  singlet  state.  In 
Section  3  we  define  correlation  functions  for  electrons  of  like  and  opposite 
spin,  and  give  the  formulas  for  these  correlation  functions  in  a  system 
of  many  electrons  in  the  AMO  approximation.  The  physical  significance 
of  these  formulas  is  discussed.  Appendix  A  contains  exact  formulas  for 
the  one  and  two  electron  density  matrices  in  terms  of  the  two  fundamental 
invariants  of  the  AMO  method.  Appendix  B  contains  a  discussion  of  the 
limited  usefulness  of  the  approximate  formulas  for  the  density  matrices 
introduced  in  Section  2. 


2.  THE  ONE  AND  TWO -ELECTRON  DENSITY  MATRICES 

The  AMO  method  involves  two  fundamental  invariants,  *  which 
are  defined  in  terms  of  two  sets  of  orthonormal  orbitals  and  a^  .  In 

a  system  of  N  =  2n  electrons,  there  are  n  orbitals  in  each  set.  The  2n 
orbitals  are  required  to  satisfy  the  further  orthogonality  relation 

^  ^KL  *  The  *unc*amental  invariants  are  p  =  anc* 
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p  =  ^KaKa*K  ‘  AMO  method  can  be  generalized  by  allowing  X 

to  depend  on  K  .  ^  In  this  case  it  is  required  that  the  2n  orbitals 
be  linearly  independent).  In  the  AMO  method  for  singlet  states,  the 
a^  are  assigned  up  spin,  the  ,  down  spin.  From  these  N 
spin  orbitals  one  constructs  the  AMO  wave  function  for  a  singlet 
state  of  an  N  electron  system  by  taking  the  singlet  projection  of  a 
determinant  of  these  spin  orbitals.  The  spinless  one  and  two  elec¬ 
tron  density  matrices  can  be  derived  from  the  general  expressions 
which  have  been  given  for  the  expectation  values  of  one  and  two  elec¬ 
tron  operators,  using  the  AMO,  singlet  state  wave  function.'*  If  we 
assume  that  \^n  «  n  *  «  1  and  neglect  as  small  all  terms  of  the 
first  order  and  higher  in  n  *  ,  we  obtain  the  following  expressions 
for  the  spinless  density  matrices; 

P,(rj')^(Ur)+^(tir)  (i) 

PzOz; W -  P,(i:r)^U;X)-f(iU)^Uir) -ft HsO^Uli) 

These  are  exactly  the  expressions  one  obtains  from  the  unprojected 
single  determinant  of  AMO's,  One  cannot  conclude  from  this  that 
for  n  -►  oo  ,  the  single  determinant  of  AMO’s  and  its  singlet  projec¬ 
tion  are  degenerate,  however.  We  discuss  this  point  in  Appendix  B. 

The  correlation  energy  in  the  Hartree-Fock  (HF)  approxima¬ 
tion  is  due  mainly  to  the  lack  of  correlation  between  electrons  of 
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opposite  spin.  Thus  we  want  to  show  that  the  AMO  method  intro¬ 
duces  such  a  correlation.  The  physical  idea  behind  the  AMO  method 
is  that  by  assigning  electrons  of  opposite  spin  to  orbitals  localized 
on  different  sublattices  of  the  alternant  system,  one  insures  that  on 
the  average  the  electrons  are  further  apart.  (An  alternant  system  is 
one  whose  lattice  can  be  constructed  from  two  equivalent,  inter¬ 
penetrating  sublattices).  In  the  AMO  method  p  is  large  about  the 
nuclei  of  one  sublattice,  p  about  the  nuclei  of  the  other  sublattice. 
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3.  ELECTRON  CORRELATION 

Before  we  can  properly  discuss  correlation  between  the  elec¬ 
trons,  we  must  discuss  the  one -electron  density  matrix.  We  have 
given  a  formula  for  the  spinless  one -electron  density  matrix.  One  might 
from  this  formula  naively  conclude  that  p  gave  the  density  of  up  spin 
electrons,  and  p  ,  the  density  of  down  spin  electrons.  This  is  just  the 
result  one  would  have  gotten  if  one  had  not  taken  the  singlet  projection 
of  the  determinant.  [We  know  of  course  that  Pj(l;  1 f )  =  Pj(  1;  t’)  + 

+  P  ^(1;  1’),  i.  e.  ,  the  one -electron  density  matrix  is  the  sum  of  an  up 

spin  density  matrix  and  a  down  spin  density  matrix.  We  have  adopted 
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the  notation  of  McWeeny  and  Mizuno  for  the  spin  components  of  the  den¬ 
sity  matrices]  .  It  is  well  known  that  in  a  singlet  state  the  spin  density 
5  +  + 

is  zero  everywhere  ,  i.e.,  Pj(is  1’)  =  Pj(l;  lf)  .  Thus  the  density  of 

up  or  down  spin  electrons  is  just  l/2(p  +  p  )  .  This  suggests  that  the 

use  of  an  unprojected  determinantal  wave  function  of  AMO's  will  give 

misleading  results  concerning  the  spin  properties  of  the  system, 
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although  it  may  give  a  quite  accurate  result  for  the  energy.  * 

For  our  purposes  it  is  sufficient  to  define  two  correlation  func¬ 
tions, 


$  (/i\ ri)  -  Pz ( *2,1 '2')  -  P,(T,f) 


%Cdj'h  -  Pz  (11;  1'Z) - p(i,i) P,  (-2 ,2) 

In  other  words,  we  say  that  two  particles  A  and  B  are  uncorrelated  when 
A  is  at  point  r^  and  B  is  at  point  r R  ,  if  the  probability  of  that  con¬ 
figuration  is  equal  to  the  product  of  the  probability  of  finding  A  at 
and  the  probability  of  finding  B  at  .  To  obtain  the  correlation  func¬ 
tions  we  have  defined,  we  need  some  spin  components  of  the  two -electron 
density  matrix.  Fortunately  for  the  singlet  state  they  are  easily  obtained 
from  the  spinless  two -electron  density  matrix  by  using  McWeeny's  and 
Mizuno’s  Eqs.  (3.2)  and  (3.5).  ^  This  leads  to  the  following  expression: 


+ 


( h;  h)=  -dz)  l(rOIOo'(z/zj-(ijll(r(i/2)lz 

<r(//z)  =p(//z)  -<5(//z) 


(2) 
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For  a  £2  »  one  sees  that  U)  r  “  ^  |  <r  ( 1 1  1)  |  <;  0  .  Thus  the 

electrons  of  opposite  spin  are  correlated  at  very  short  distances.  One 
feels  that  this  is  the  energetically  important  correlation.  From  the 
unprojected  determinant  one  finds  1 2;  12)  sO  ,  i.e. ,  the  electrons 

of  opposite  spin  are  uncorrelated  in  a  formal  sense,  although  the 
energy  per  electron  in  a  many  electron  system  is  approximately  equal 
to  that  of  the  singlet  projection.  One  can  argue  that  to  discuss  correla¬ 
tion  between  electrons  of  opposite  spin  using  an  unprojected  determinantal 
wave  function  is  not  meaningful  since  the  eigenfunctions  of  a  spinless 
Hamiltonian  must  be  spin  eigenfunctions.  This  is  consistent  with  our 
observation  on  the  properties  of  the  one -electron  density  matrix. 

In  addition  to  the  short  range  correlation  between  the  electrons, 
there  is  also  a  long  range  correlation.  For  |r^  *  £3 1  80  large  that 

p(l|2)  and  p(  1 1 2)  are  negligibly  small  in  comparison  to  p(  1 1  l)  and 
p(  1 1  i),  the  correlation  functions  are  related, 

Cjz  ( 12,  ]~2)  s?  -(/2)'<7(//i)  <1(2.1 2)  =  -9J/2;  t2)  (3) 

The  physical  meaning  of  this  equation  is  easily  extracted  by  recalling  the 
properties  of  p  and  p  in  alternant  systems.  If  we  choose  r^  so  that 
p(  1 1  1 )  is  a  maximum,  g^(l2;tZ)  has  a  maximum  wherever  p( Z |  2) 

has  a  maximum,  and  (lx;  i2)  has  a  maximum  wherever  p(Z  |  Z)  has 
a  maximum.  If  r.  is  chosen  so  that  p(l|  l)  is  a  maximum,  then 

+  -  ~~l  .  44  44 

g-(l2;  12)  has  a  maximum  wherever  p(2|2)  haB  a  maximum,  and  g~(l2;  12), 
wherever  p(Z  ]  2)  has  a  maximum.  Thus  in  this  relative  sense,  electrons 
of  one  spin  are  localized  on  one  sublattice  of  an  alternant  system,  electrons 
of  the  other  spin  on  the  othex*  sublattice. 

Although  we  have  explicitly  used  the  AMO  approximation  in  the 
preceding  discussion,  similar  conclusions  would  be  reached  if  one  took 

g 

the  singlet  projection  of  an  unrestricted  Hartree-Fock  (UHF)  wave  function. 
In  that  case,  however,  one  would  not  be  able  to  say  in  general  that  tie 
electrons  of  opposite  spin  were  localized  on  different  sublattices  of  the 
system,  i.e.,  you  would  need  information  about  the  p  and  p  of  the  UHF 
approximation . 

In  conclusion,  we  remark  that  the  results  of  this  paper  provide  a 

physical  understanding  of  the  AMO  wave  functions  which  the  author  has  coa- 

9 

sidered  for  the  electron  gas. 
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3.  ELECTRON  CORRELATION 


Before  we  can  properly  discuss  correlation  between  the  elec¬ 
trons,  we  must  discus*  the  one -electron  density  matrix.  We  have 
given  a  formula  for  the  spinless  one -electron  density  matrix.  One  might 
from  this  formula  naively  conclude  that  p  gave  the  density  of  up  spin 
electrons,  and  p  ,  the  density  of  down  spin  electrons.  This  is  just  the 
result  one  would  have  gotten  if  one  had  not  taken  the  singlet  projection 
of  the  determinant.  [We  know  of  course  that  Pj(l;  l1)  =  P^(  1;  1')  + 

+  P  v ( I ;  l1),  i.e.,  the  one -electron  density  matrix  is  the  sum  of  an  up 
spin  density  matrix  and  a  down  spin  density  matrix.  We  have  adopted 

5 

the  notation  of  McWeeny  and  Mizuno  for  the  spin  components  of  the  den¬ 
sity  matrices] .  It  is  well  known  that  in  a  singlet  state  the  spin  density 
is  zero  everywhere  ,  i.e.  ,  Pj(i;  1')  =  Pj(I»  l1)  •  Thus  the  density  of 
up  or  down  spin  electrons  is  just  l/2(p  I  p)  .  This  suggests  that  the 
use  of  an  unprojected  determinantal  wave  function  of  AMO's  will  give 

misleading  results  concerning  the  spin  properties  of  the  system, 

6  7 

although  it  may  give  a  quite  accurate  result  for  the  energy.  * 

For  our  purposes  it  is  sufficient  to  define  two  correlation  func¬ 
tions, 


-h  ’h 


PJlZ-,1'2')-  P,(l,l)  P,U;2‘) 


-  Pjn^-Pd'hpu^) 


+  + 


+ 


In  other  words,  we  say  that  two  particles  A  and  B  are  uncorrelated  when 
A  is  at  point  r^  and  B  is  at  point  r ^  ,  if  the  probability  of  that  con¬ 
figuration  is  equal  to  the  product  of  the  probability  of  finding  A  at 
and  the  probability  of  finding  B  at  ,  To  obtain  the  correlation  func¬ 
tions  we  have  defined,  we  need  some  spin  components  of  the  two -electron 
density  matrix.  Fortunately  for  the  singlet  state  they  are  easily  obtained 
from  the  spinless  two-electron  density  matrix  by  using  McWeeny's  and 
Mizuno’s  Eqs.  (3.2)  and  (3.5).  ^  This  leads  to  the  following  expression; 


( 11;  7z)~  -  ( lzf(T(lh) <r(z/2j -(C) 'l 0~(//jz)IZ 

<r(iU)  1 1 z)  ~fOUJ 


(2) 
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F°r  £|  =  £2  *  one  sees  that  g2(ll;  ll)  =  -  ^  | <r(l  |  l)| ^  <  0  .  Thus  the 

electrons  of  opposite  spin  are  correlated  at  very  short  distances.  One 

feels  that  this  is  the  energetically  important  correlation.  From  the 

♦  + 

unprojected  determinant  one  finds  g 1 2 ;  l2)  s  0  ,  i.e.  ,  the  electrons 
of  opposite  spin  are  uncorrelated  in  a  formal  sense ,  although  the 
energy  per  electron  in  a  many  electron  system  is  approximately  equal 
to  that  of  the  singlet  projection.  One  can  argue  that  to  discuss  correla¬ 
tion  between  electrons  of  opposite  spin  using  an  unprojected  determinantal 
wave  function  is  not  meaningful  since  the  eigenfunctions  of  a  spinless 
Hamiltonian  must  be  spin  eigenfunctions.  This  is  consistent  with  our 
observation  on  the  properties  of  the  one -electron  density  matrix. 

In  addition  to  the  short  range  correlation  between  the  electrons, 
there  is  also  a  long  range  correlation.  For  (r^  -  r^l  so  large  that 
p(l|2)  and  p(l|2)  are  negligibly  small  in  comparison  to  p(l|  l)  and 
p(l  |  l),  the  correlation  functions  are  related# 

( 12,  U)  =  -(llj'vf//!)  <T(2.lz)  =  -gjtei  ffl  (3) 


The  physical  meaning  of  this  equation  is  easily  extracted  by  recalling  the 
properties  of  p  and  p  in  alternant  systems.  If  we  choose  r  ^  so  that 
p(l|l)  is  a  maximum,  g~(12;  t2)  has  a  maximum  wherever  p(2|2) 
has  a  maximum,  and  (12;  12)  has  a  maximum  wherever  p (Z  |  E)  has 
a  maximum.  If  is  chosen  so  that  p(l|  l)  is  a  maximum,  then 
g2(l2;  12)  has  a  maximum  wherever  p(2|2)  has  a  maximum,  and  g 2(  1 2 ;  12), 
wherever  p(2|2)  has  a  maximum.  Thus  in  this  relative  sense#  electrons 
of  one  spin  are  localized  on  one  sublattice  of  an  alternant  system,  electrons 
of  the  other  spin  on  the  other  sublattice. 


Although  we  have  explicitly  used  the  AMO  approximation  in  the 
preceding  discussion,  similar  conclusions  would  be  reached  if  one  took 

8 

the  singlet  projection  of  an  unrestricted  Hartree-Fock  (UHF)  wave  function. 
In  that  case,  however,  one  would  not  be  able  to  say  in  general  that  the 
electrons  of  opposite  spin  were  localized  on  different  sublattices  of  the 
system,  i.e.,  you  would  need  information  about  the  p  and  p  of  the  UHF 
approximation . 

In  conclusion,  we  remark  that  the  results  of  this  paper  provide  a 

physical  understanding  of  the  AMO  wave  functions  which  the  author  has  con- 

o 

sidered  for  the  electron  gas. 
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APPENDIX  A 

Exact  Formulas  for  P^  and  P^  for  the  Singlet  State 
In  addition  to  p  and  p  ,  we  introduce  the  functions 


CO  (HZ)  -  fdr2  (C  (U  3)  f  (3l2)}  (a  i) 

LO(l/2)  »  [_Uj(Zl l)J*  (A2) 


The  spinless  one -electron  density  matrix  is 

P,n,n  =f,w)[? dliOtf  (III1)  i-  fjjtf W//, !•)-<■& or t')J  (A3) 
f,(A2)  -  n'd-P)'1  (i-T'T In  -(n+i)Y+X”z3  m 

ft(7f )  -  n'(i-?J'(i-r“Y,ll-(nOX'+n7r1]  m 


The  spinless  two-electron  density  matrix  is 

P//2,/2'J-0V{[ X' I'd H ('C +? (212)]  - 


+[?  oil) +f(t, /  nj[ 0212'}  *  u>(2/2')J 

-t [cj(i//)-f(Z)(///l)][^(2l2'J  *  q(z/z')] 

-  I<z(il2!)+p(il2!)][(jz/i')  +od(2/rj] 

-  [uX  111 )  t  U)(ll2!)][q(z/l'J  i-  f  (2! / ')] 

•+  00(1121)  CO  (zll  ')  +  U)  (liz  ')oj(Z  1 1') J 

-+  fs(^z){  (jh(Hi)lo(zIzI) + uo( I ! i)uj(zlz') 

-  iuXl/2')u)(2/l)  CO(//z')cD(Z/r)J 
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fjy)  -  n'u-tru-n1  d-^y'x. 

x  [n(n-i) ~2(n-l)X~i n(tni)'Kl-2.  X n'2J 
(jy>  ~  ri'(n-/)-‘(i~tf)~2( i-Azn^'r>  x  ( 

*f  »-/  ~(rmn*-+  (n-fOXn-(n-l)r^J 

fs(AV -  4n'(n-ir‘( I- i 
x  [2 - n(rui)'X'-z+2(r?--l)'Xn-n(r,-i)Xn~] 


APPENDIX  B 

Validity  of  Formulas  for  P^  and  for  the  Case  of  Large  n 
Assuming  thai  X^n  «  n  *  «  1  ,  one  finds  to  order  n  ^  , 

-£(>.2)-  i-rr'd-XY'X^ 

{^‘n'd-W1  (B1) 

40M  -  l-Zn-'U-XT"f+ihHi-X)~zX' 

U  Qf) = rtd~  ^r-  i  n-z(/- rr2  a  * 

-  4n~z(i -)?)'*■  <B2> 


- 1  -2 

If  we  write  Pj  as  a  power  series  in  n  ,  then  to  order  n 


P,  =  ^'0J+ 


-/ 


(B3) 


(In  Eq.  (1)  we  took  the  limit  n  -►  00  ,  so  that  the  P^  given  there  corre¬ 
sponds  to  P  ^  in  the  above  equation.)  From  Eq.  (A3)  we  see  that  term 
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by  term  p/  i a  proportional  to  the  electron  density.  Integrated  over 

1  (1) 

all  space  P ^  9  gives  zero,  as  is  implied  by  the  expansion  in  n 
(The  trace  of  is  equal  to  2 n  ).  Thus  corresponds  to  a 

charge  distribution  of  zero  net  charge.  To  the  expectation  value  of  the 
nuclear  coulomb  potential,  n  gives  a  contribution  independent 

of  n  since  the  number  of  nuclei  is  proportional  to  n  .  (If  the  trace  of 
Pj^  were  proportional  to  n  as  a  term  by  term  inspection  of  P^^ 
would  imply,  then  it  would  contribute  to  the  expectation  value  of  the  sum 
of  the  nuclear  coulomb  potentials,  a  term  proportional  to  n  .  On  the 

other  hand  to  the  trace  of  P,  it  would  contribute  a  term  independent  of 

1  (0) 
n  ,  which  would  be  negligibly  small  in  comparison  to  the  trace  of  Pv  9 

for  large  n  .  This  suggests  that  the  direct,  numerical  calculation  of 

density  matrices  for  many  electron  systems  could  give  expectation  values 

in  error  by  terms  proportional  to  the  total  number  of  electrons,  although 

the  error  in  normalization  would  be  negligibly  small). 

A  similar  analysis  may  be  carried  through  for  P7  .  In  that  case 

(l)  L 

the  part  of  P^v  9  which  is  term  by  term  of  the  order  of  magnitude  of  the 

electron  density  squared,  can  be  factored  into  two  parts.  The  one  part 

is  a  charge  distribution  of  zero  net  charge,  the  other  has  a  net  charge 

proportional  to  n  .  Thus  P^*^  contributes  to  the  total  energy  a  term 

independent  of  n  for  the  singlet  state. 
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FOOTNOTES 
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5.  R.  McWeeny  and  Y.  Mizuno,  Proc.  Roy.  Soc.  (London)  A 2_59,  554 
(1961). 

6.  Criticisms  of  the  use  of  the  unrestricted  Hartree-Fock  method  in 
computing  spin  densities  have  been  given  recently  by  R.F.  Wood  and 

K. F.  Berggren  (Phys.  Rev.  to  be  published). 

7.  See,  however,  the  discus sion  given  by  W -  Marshall,  Proc.  Phys. 
Soc.  7jS,  113  (1961). 

8.  F.  Sasaki  (private  communication)  and  P.O.  Ldwdin  (private 
communication) . 

9.  One  such  wave  function  has  been  described  briefly  in  Section  5  of  a 
paper  by  the  author.  [Phys.  Rev.  127 ,  1650  (1962)  ]. 
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